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Abstract 

We consider a general model for representing and manipulating parametric curves, in which 
a curve is specified by a black box mapping a parameter value between and f to a point in 
Euclidean d-space. In this model, we consider the nearest-point- on- curve and farthest-point- on- 
curve problems: given a curve C and a point p, find a point on C nearest to p or farthest from p. 
In the general black-box model, no algorithm can solve these problems. Assuming a known 
bound on the speed of the curve (a Lipschitz condition), the answer can be estimated up to an ad- 
ditive error of e using 0{1/ e) samples, and this bound is tight in the worst case. However, many 
instances can be solved with substantially fewer samples, and we give algorithms that adapt to 
the inherent difficulty of the particular instance, up to a logarithmic factor. More precisely, if 
OPT(C, p, e) is the minimum number of samples of C that every correct algorithm must perform 
to achieve tolerance £, then our algorithm performs 0(OPT(C,p, e) log(e~^/OPT(C,p, e))) sam- 
ples. Furthermore, any algorithm requires 0(fc log(e~^/fc)) samples for some instance C" with 
OPT(C",p, e) = k] except that, for the nearest-point-on-curvc problem when the distance be- 
tween C and p is less than e, OPT is 1 but the upper and lower bounds on the number of 
samples are both Q{l/£). When bounds on relative error are desired, we give algorithms that 
perform 0(OPT • log(2 + (1 + e^^) • jti^^/OPT)) samples (where m is the exact minimum or 
maximum distance from p to C) and prove that ri(OPT • log(l/e)) samples are necessary on 
some problem instances. 



1 Introduction 

Computational geometry has traditionally been focused on polygonal objects made up of straight 
line segments. In contrast, applications of geometric algorithms to computer-aided design and 
computer graphics usually involve more complex curves and surfaces. In recent years, this gap 
has received growing attention with algorithms for manipulating more general curves and surfaces, 
such as circular arcs |DFMTn2] . conic arcs |BEH+n2irWei?i2j . and quadratic surfaces The 
most general type of curve commonly considered in this algorithmic body of work is a piecewise 
bounded-degree polynomial (algebraic) curve, although such curves are not usually manipulated 
directly and are more typically assumed to govern some process such as the motion of a polygon 



in kinetic collision detection |BEG"'"99 



Parametric Black-Box Curves. A much more general model for specifying curves is the para- 
metric black-box model that represents a curve in Euclidean d-space as a function C: [0, 1] — > M'^. 
The only operation that can be performed is to sample (evaluate) the function at a given parameter 
value X £ [0, 1]. 
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Figure 1: An instance of the nearest-point-on-curve problem. 



Solving any nontrivial problem on a black-box curve requires some additional conditions on the 
behavior of the curve. We assume the Lipschitz condition that ||C(xi) — C(x2)|| < L\xi — X2\ for 
all xi,X2 G [0,1], for a known constant L. Any piecewise-C^ curve has such a parameterization. 
By uniformly scaling the curve in W^, we can assume that the Lipschitz constant L is 1. 

Nearest- and Farthest-Point-on-Curve Problems. In this paper, we solve two of the most 
basic proximity queries about black-box Lipschitz curves: given a curve C and a point p, find 
a point on C that is closest to p {nearest point), and find a point on C that is farthest from p 
(farthest point). In the black-box model, these problems are impossible to solve exactly, because an 
algorithm will never, in general, sample the nearest or farthest point. Thus, a problem instance also 
specifies an additive error tolerance e, and our goal is to find a point on the curve C whose distance 
to the point p is within ±e of the minimum or maximum possible. See Figure ID Although we focus 
on absolute (additive) error in this paper, we show in Section El how to modify the absolute-error 
algorithms to obtain relative-error algorithms (whose output is accurate to within a factor of e) 
that have nearly optimal adaptive performance. 

Hard and Easy Instances. Any nearest-point-on-curve or farthest-point-on-curve instance can be 
solved using l/2e + 0(1) samples: C(0), C(2e), C(4e), . . . , C(l). Unfortunately, this many samples 
can be necessary in the worst case. For example, when C{x) = q for all x outside an interval of 
length 2e where at speed 1 the curve moves toward p and then returns to q, we need l/2e — 0(1) 
samples to find the interval. Thus, worst-case analysis is not very enlightening for this problem. 

On the other hand, many instances are substantially easier. As an extreme example, if C is a 
unit-length line segment, then two samples, at C(0) and C(l), completely determine the curve by 
the Lipschitz condition. 

Adaptive Analysis. Because the instance-specific optimal number of samples varies widely from 
0(1) to 0(l/e), we use the adaptive analysis framework, considered before in the context of boolean 
set operations |DL()Mfln] as well as sorting |E(]W92j and aggregate ranking |FLNn3j . In the adap- 
tive analysis framework, the performance of an algorithm on a problem instance is compared to 
OPT, the performance of the best possible algorithm for that specific problem instance. By defini- 
tion, for every problem instance, there exists an algorithm that achieves OPT on that instance. The 
question is whether one adaptive algorithm uses roughly OPT(C, p, e) samples for every instance 
iC,p,e). 
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Our Results. We develop adaptive algorithms that solve the nearest-point-on-curve and farthest- 
point-on-curve problems using 0{OFT{C,p,£) log(e-VOPT(C,p,e))) samples; except that, for the 
nearest-point-on-curve problem when the distance between C and p is less than e, the number of 
samples may be Q{l/e), yet OPT = 1. We also prove that these algorithms are optimally adaptive 
in the sense that no adaptive algorithm can achieve a strictly better bound (up to constant factors) 
with respect to OPT and e. Specifically, we show that, for any e > and k > 0, there is a family 
of curves C each with OPT(C, p, e) = k such that every algorithm (even randomized) requires 
^{klog{e~^ /k)) samples on average for a curve C selected uniformly from the family; and there is 
a family of instances of the nearest-point-on-curve problem where the distance between C and p is 
less than e such that every algorithm requires ^}{l/e) samples on average, but OPT is 1. 

Related Work. Because our curve model is a black box, the problems that we consider here 
have natural formulations in information-based complexity terms (see |TWW88| for an overview). 
However, information-based complexity is primarily concerned with worst-case, average-case, or 
randomized analysis of more difficult problems, rather than adaptive analysis of algorithms for 
easier problems as in this paper. Information-based complexity does consider adaptive algorithms 
(algorithms for which a query may depend on the answers to previous queries) , but primarily when 
they are more powerful than non-adaptive algorithms in the worst (or average, etc.) case, such as 
for binary search. 

The problem of maximizing a Lipschitz function has been studied in the context of global 
optimization. This problem essentially corresponds to the special case of the nearest- or farthest- 
point-on-curve problem in which d = 1. Beyond worst-case analysis, many algorithms for this 
problem have been studied only experimentally (see, e.g. [HJ95]), but Piyavskii's algorithm |Piy72| 
has been previously analyzed in what is essentially the adaptive framework, first in |Dan71j . The 
analysis was sharpened in |HJL91j to show that the number of samples the algorithm performs on 
(C, e) is at most 4 times OPT(C, e). As Theorem El shows, this analysis cannot generalize to d > 1. 

Practitioners who manipulate curves and surfaces typically use numerical algorithms, which 
are extremely general but sometimes fail or perform poorly, or specialized algorithms for specific 
types of curves, such as B-splines. Some algorithms for manipulating general parametric curves and 
surfaces guarantee correctness, but the theoretical performance of these algorithms is either not 
analyzed |Sny92 or analyzed only in the worst-case [jTiQS , GW90^ . At the heart of our algorithm is 



Giinther and Wong's jOWDOj observation that the portion of a Lipschitz curve between two nearby 
sample points can be bounded by a small ellipse, as described in Section f3. II 



2 Problem Statement 

We use the real RAM model, which can store and manipulate exact real numbers in 0(1) time 
and space. Manipulation of real numbers includes basic arithmetic (+, — , x, -^), comparisons, and 
nth roots. We separately analyze the number of samples and the additional computation time. 
Although we describe and analyze our algorithms in R^, both the algorithms and their analyses 
trivially carry over to for d > 2. 

We assume without loss of generality that the Lipschitz constant is 1 , that the parameter space 
is the unit interval, and that the query point p is the origin O. Throughout our discussion of 
nearest-point-on-curve, dmin refers to the minimum distance from C to the origin. Analogously, 
dmax denotes the maximum distance from C to the origin. We assume that e is smaller than 1/2 
because otherwise, a single sample at 1/2 immediately solves both problems. The two problems we 
consider are to find a point on C whose distance to O is approximately dmin or dmax- 
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Figure 2: Some possible curves C inside an ellipse. 



Problem Nearest-Point-On-Curve Given a Lipschitz curve C and an < e < 1/2, find a 
parameter x such that ||C(2;)|| < dmin + 

Problem Farthest-Point-On-Curve Given a Lipschitz curve C and an < e < 1/2, find a 
parameter x such that ||C(2;)|| > (imax — 



3 Nearest-Point-On-Curve: Adaptive Algorithm and its Analysis 
3.1 Main Idea 

The main observation is that, if we have sampled C at xi and X2, then for x between xi and X2, 
||C(x) — C(xi)|| + ||C(x) — C(x2)|| < |x2 — xi|. This means that when the parameter x is between 
xi and X2, C{x) stays within an ellipse with foci at C(xi) and C(x2), whose major axis (sum of 
distances to foci from a boundary point) has length |x2 — xi|. See Figure |21 Note that this ellipse 
is tight: by changing C only between xi and X2, we can force it to pass through any point in the 
ellipse while keeping C Lipschitz. The following propositions formalize this idea: 

Definition 1. Given a Lipschitz curve C and an interval [xi,X2] ^ [0, 1], define the ellipse 

Ec{xi,X2) = [p (^^^ I ||C7(xi)-p|| + ||C7(x2)-p|| <X2-Xi}. 
Proposition 1. For an interval J = [xi,X2] ^ [0, 1], C(J) C Ec{xi,X2)- 

Proof: Let x G J. By the Lipschitz condition, ||C(xi) — C(x)|| < x — xi and similarly ||C(x2) — 
C(x)|| < X2 - X. Adding these, we get ||C(xi) - C(x)|| + ||C(x2) - C(x)|| < X2 - xi, so C(x) G 
Ec{xi,X2). □ 

Proposition 2. Let J = (xi,X2) ^ [0, 1] and let C he a Lipschitz curve. Then for every point p 
in Ec{xi,X2), there is a Lipschitz curve C such that C(x) = C"(x) for x J and for some x £ J, 
C'{x) =p. 

Proof: We can make C on J consist of a line segment from C(xi) to p and another one from p to 
C(x2). Because the total length of these line segments is at most X2 — xi, we can parametrize C 
at unit speed (or less) on J. □ 

The following proposition will often be used implicitly in our reasoning: 

Proposition 3. If J' = [x'j^jXg] and J = [xi,X2] and J' C J, then Ec{x'i,X2) ^ Ec{xi,X2)- 
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Proof: If p G Ec{x[,X2) then ||C(2;'^) — p\\ + ||C(x2) — p\\ < X2 — x[ by definition. We have 
||C(xi) — C(x'^)|| < x'l — xi and ||C(x2) — C(x2)|| < X2 — x'2 by the Lipschitz condition on C. 
Adding the three inequahties and applying the triangle inequality twice, we get: 

\\C{x^)-v\\^\\v-C(x2)\\ < 

< ||c(xi) - c{x[)\\ + \\c{x[) - p\\ + Hp - + ||C(4) - cix2)\\ <X2- x,. 

So p £ Ec{xi, X2), as required. □ 

For notational convenience, let closest-possible{xi, X2) denote the minimum distance from a 
point in Ec{xi, X2) to the origin. 

3.2 Proof Sets 

The properties of Eq immediately suggest a criterion for determining whether a set of points on 
a curve is sufficient to guarantee that a point sufficiently close to O is among those in the set: 
the distance from O to the nearest sampled point and the distance from O to the nearest ellipse 
(around adjacent points) should differ by at most e. 

Definition 2. Let P = {xi, X2, ■ ■ ■ , Xn} be a set of parameters in [0,1] so that = xi < X2 < 
• • • < Xn = 1. Let Xmin € P be an element that minimizes ||C(xi)||. Then P is a proof set if 
\\C{xmin)\\ — £ < closest-possible(3;j, Xj+i) for all i. 

The following proposition shows that producing a proof set is the only way an algorithm can 
guarantee correctness. 

Proposition 4. Let P = {xi, X2, . . . , x„} C [0, 1] so that = xi < X2 < ■ ■ ■ < x„ = 1. Let x^iin 
be an element of P that minimizes ||C(2;j)||. If P is a proof set, then for any curve C such that 
Xrmn is a Solution to nearest-point- on- curve. Conversely, if P is not a proof set, 
there is a curve C such that C'{xi) = C{xi) for all i and for which Xmin is not a solution. 

Proof: For any curve C for which C'{xi) = C{xi), P is a proof set for C' precisely when it is a 
proof set for C. Applying Proposition ^ we find that C([0, 1]) is contained in the union of the 
ellipses Ec{xi,Xi+i). So, if P is a proof set, ||C(xmin)|| — £ < l|C"(x)|| for all x G [0,1], which 
implies that Xmin is a solution for C 

Conversely, if P is not a proof set, then there is a point p in some ellipse £'c(xj, Xj+i) such that 
||C'(2;min)|| — s > \\p\\- By Proposition 12 we can construct a curve that coincides with C except in 
and passes through p. For this curve, Xmin will not be a solution. □ 

The requirement that xi = and x„ = 1 allows the analysis to avoid special cases. An 
algorithm could guarantee correctness without sampling these endpoints, but because this saves 
only a constant amount of work, we ignore this possibility in favor of simpler analysis. 

3.3 Algorithm Description and Correctness 

As we sample the curve, we maintain a set of ellipses around the unsampled intervals. At each 
step, we take the interval whose ellipse is closest to the origin and sample in the middle of it, thus 
replacing it with two smaller intervals (with smaller ellipses). When the sampled points form a 
proof set, we terminate and output the closest point of those sampled. 
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Let Q be a priority queue that stores triples of real numbers {d,xi,X2) sorted by d. The 
algorithm is as follows: 

closest-point(C, e) 

1. Add (closest-possible(0, 1), 0, 1) to Q 

2. If ||C(0)|| < ||C(1)|| then (x„,in,dmin) ^ (0, ||C(0)||) 
else (Xmin,C^mm) ^ (1, || C(l) || ) 

3. Do until finished: 

4. {d,xi,X2) ^ extract-min((5) 

5. If dmin — £ < d then output Xmin and STOP 

6. X ^ (xi + X2)/2 

7. If ||C(x)|| < drain then (Xmin,dmin) ^ (x, ||C(x)||) 

8. Add (closest-possible( ) to Q 

9. Add (closest-possible(x, 3:2), X, X2) to Q 

Correctness follows from Proposition \^ the algorithm stops when the points sampled form a 
proof set and outputs the closest point. To show termination, we note that no interval of length 2e 
or less is ever subdivided: 

Proposition 5. // in line 5, X2 — xi < 2e, closest-point terminates at this line. 

Proof: Because dmm stores the minimum known distance to a point, dmin < l|C'(xi)|| and dmm < 
||C(x2)||. Let p be a point in Ec{xi,X2) such that ||p|| = d. Then by the definition of Ec, 
\\C{xi) —p\\ + ||C(x2) —p\\ < 2e. This means that at least one of ||C(xi) —p\\ < e or ||C(x2) —p\\ < e 
is true. If ||C(xi) — p\\ < £, then, by the triangle inequality, ||C(3;i)|| — \\p\\ < e. This implies that 
dmin — d < £ so the algorithm stops. Similarly for the other possibility. □ 

Prom this proposition, we can conclude that closest-point stops after at most 0(l/e) loop 
iterations because only 0(l/e) sample points at least e apart can fit in [0, 1], and in each iteration 
of the loop, the algorithm always samples one new point in an interval of width at least 2e. 

3.4 Ellipse Lemma 

To analyze CLOSEST-POINT, we will make use of one geometric fact in three incarnations: 

Ellipse Lemma. Let < xi < X2 < X3 < X4 < 1. Also, let d, a G M with < a < d. 
// closest-possible(xi, X2) < d, closest-possible(x3, X4) < d, and ||C(x2)|| > d + a, then closest- 
possible(xi, X4) < d — a. 

Proof: See PigurelHl We may assume without loss of generality that X2 = X3, because if closest- 
possible(x3, X4) < d, then closest-possible(x2, X4) < d. Let p be the intersection of the circle 
\\v\\ = d — a and the ray from the origin through C(x2). Obviously, ||p|| = d — a. We will show 
that ||C(xi) — p\\ < X2 — xi and ||C(x4) — p\\ < X4 — X2. This will prove that p £ Ec{xi,X4), and 
therefore closest-possible(xi, X4) < d — a. 

Because closest-possible(xi, X2) < d, there is a point q such that ||g|| < d and ||C(xi) — f?!! + — 
C(2;2)|| < X2 — xi. We may set the axes so that C(x2) is on the y axis. So let C(x2) = (0, y). Then 
y > d+a and p = {0,d — a). Now if g = {xq,yq), then — C(x2)|| = y^'x^ + {y — yqY and Hg— p|| = 

\Jxl + (yq - {d- a))2. Because \\q\\ < d, yg < d, we have {yq-{d-a)f < {{d+a)-yqf < {y-yqf, 
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Figure 3: Ellipse Lemma 



which means that H? — p|| < \\q ~ C{x2)\\. Using this, the triangle inequality, and the construction 
requirement of g, we get: 

\\C{xi)-p\\ < \\C{xi)-q\\ + \\q-p\\ < 

< \\C{xi) - q\\ + \\q - C{X2)\\ <X2- XI. 

The argument that ||C(x4) — p\\ < 2:4 — X2 is symmetric. □ 

When generalizing this lemma from to W^, we consider separately the planes through 
O, C{xi),C{x2) and through O, C(x2), C(x4). 

We will use the Ellipse Lemma in three different places in the analysis, so we prove three simple 
corollaries: 

Ellipse Lemma (1). Let [xi,x] and [x,X2\ he intervals. LetO < e < d. /f closest-possible(3;i, X2) > 
d — e and ||C(2;)|| > d, i/ien closest-possible(xi, x) > d — e/2 or closest-possible(x, a;2) > d — e/2 or 
both. 

Proof: Suppose for contradiction that this is not the case: that we have closest-possible(xi, x) < 
d — e/2 and closest-possible(x, X2) < d — e/2. Let d' be the larger of closest-possible(xi, x) and 
closest-possible (x, X2). We have d' < d — e/2. Now let a = d — d' and apply the Ellipse Lemma to 
[xi, x], [x, X2], d' and a to get closest-possible(xi, X2) < d' — a. But d' — a = 2d' — d < d — e, which 
contradicts the assumption that closest-possible(xi, X2) > d — e. □ 
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Ellipse Lemma (2). Let [xi,X2] and [x3,X4] be intervals with > X2- Let < e/2 < d. If 
closest-possible(a;i, 3:2) < d, closest-possible(a;3, X4) < d, and ||C(x2)|| > d + e/2, then closest- 
possible(xi, X4) < d — e/2. 

Proof: Let a = ||C(x2)|| — d so a > e/2. Now apply the Ellipse Lemma to [xi,a;2], [x3,X4], d, and 
a to get closest-possible(rEi, 2:4) < d — a < d — e/2. □ 

Ellipse Lemma (3). Let [xi,x] and [x,X2] be intervals and suppose < e/2 < d. /f closest- 
possible(xi, X2) > d — e/2 and ||C(x)|| > d + e/2, then closest-possible(xi, x) > d or closest- 
possible(x, X2) > d or both. 

Proof: Suppose that closest-possible(xi, x) < d and closest-possible(x, X2) < d. Apply the Ellipse 
Lemma(2) to get that closest-possible(xi, X2) < d — e/2, which is a contradiction. □ 

3.5 OPT 

We define the OPT of a problem instance to be the number of samples that the best possible 
algorithm makes on that instance. In our analysis, we use the fact that OPT is equal to the 
size of the smallest proof set, which follows from Proposition^ Note that OPT depends on C 
and e, but we write OPT or OPT(e) instead of OPT(C, e) when the arguments are clear. For 
the analysis of closest-point with e < dmin we need the following estimate: for any curve C, 
0PT(e/2) = 0(OPT(e)). We prove this by starting with a proof set for e, inserting a new sample 
point in between every pair of sample points in the proof set, and using the Ellipse Lemma with a 
continuity/connectedness argument to show that we can force the result to be a proof set for e/2. 

Proposition 6. If e < dmin. for any problem instance (C,e), 0PT(C,e/2) < 20PT(C,e). 

Proof: Consider a proof set P of size OPT(e). Let Xj be the i*^ smallest element of P. Because 
P is a proof set, closest-possible(xj, Xj+i) > dmin — Now let 

A = {x £ [xj,Xj+i] I closest-possible(xj, x) > dmm — s/2} 
B = {x £ [xj,Xj-|_i] I closest-possible(x, Xj+i) > dmm — s/2} 

By the Ellipse Lemma(l), AUB = [xj,Xj+i]. Also, because closest-possible(xj, Xj) = ||xj|| > dmin 
and closest-possible(xi+i, Xj+i) = ||xj+i|| > dmm, A ^ % and i? / 0. Because closest-possible is 
continuous in both variables, A is closed relative to [xi,Xi+i], being the preimage of the closed set 
{t \ t > dmin — e} under t = closest-possible(xj, x) with respect to the second variable. Similarly, 
B is closed relative to [xi,Xi+i]. Because [xj,Xj+i] is connected, A n / 0, so let x S j4 n i?. 
This means closest-possible(xj, x) > dmin — e/2 and closest-possible(x, Xj+i) > dmin — e/2. So for 
every pair of adjacent samples in P, we can insert a new sample x between them (x may, of course, 
coincide with one of the samples already in P, in which case we ignore it) so that in the resulting set, 
closest-possible(xj, Xj+i) > dmin — e/2 for all j. Thus, we will have inserted at most an additional 
|P| — 1 elements. In order to make the result a proof set, we may need to insert one more element, 
X such that ||C(x)|| = dmin- This will make the result into a proof set for e/2 with 2|P| elements. □ 

3.6 Phases 

We split an execution of closest-point into two phases and analyze each phase separately, giving 
an upper bound on the number of curve samples. The phases are a construction for the analysis 
only; the algorithm does not know which phase it is in. The algorithm starts out in Phase 1, and 
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switches to Phase 2 when all of the ellipses around intervals stored in Q are no closer than dmin — e/2 
to the origin. The distance from the ellipses in Q to the origin can only grow (as ellipses close to 
the origin are replaced by ellipses farther away), so once the algorithm enters Phase 2, it can never 
leave it. Let P be a proof set for e/2 whose size is 0PT(e/2). We show that in each phase, the 
number of samples is 0(|P| log(e~^/|P|)). We will want the following easy fact: 

Proposition 7. Let a-i for 1 < ^ < |P| be positive real numbers. If we have X^^^^ Oi < e~ , then 
EE'iloga, < |P|log(e-V|P|). 

P| / ipi 5^-^!' ai -1 

Proof: By the arithmetic-geometric mean inequality, we have y JJ-^^^ flj < — < jp|-. Taking 

the logarithm of both sides gives us |p°^ ^ log(e~^/|P|). Multiplying both sides by |P| gives 
us the desired result. □ 



3.7 Phase 1 

Proposition 8. // closest-possible(xi, X2) < c^min — e/2, then P must have a point in the open 
interval (xi, X2). 

Proof: For contradiction, suppose that P has no point in (xi,X2)- This means that P has two 
consecutive points, x'^ and X2, such that [xi,X2] C [x'^,a;2]. So Ec{xi,X2) C Ec{x'i,X2) and there- 
fore, closest-possible(xi, X2) > closest-possible(x']^, X2), which means that closest-possible(x']^, 0:2) < 
dmin — e/2. Hence, P cannot be a proof set for e/2. □ 

Let J = [xi, X2] ^ [0, 1] be an interval that is subdivided in Phase 1. This implies that closest- 
possible(xi, X2) < c^min — e/2 so by Proposition |H1 P must have a point in {xi,X2). This means that 
any interval that is subdivided in Phase 1 contains a point of P. 

We need to count the samples in this phase. We achieve this by classifying every subdivision 
as either a "split" or a "squeeze" . A subdivision is a split if both resulting intervals contain points 
from P and a squeeze if one of the resulting intervals has no points from P. Because the number 
of splits cannot be more than |P| — 1, we only need to count squeezes. If J is an interval in Q at 
some point in the execution of Phase 1, let S{J) be the number of squeezes that have happened to 
intervals containing J and let L{J) be the length of J. We want the following invariant: 

Proposition 9. If at some point during Phase 1 of the algorithm, the intervals that intersect P 
are Ji,J2,...,Jk, then Zti '^^^^'^L{Ji) = 1. 

Proof: We proceed by induction on the number of subdivisions. At the start of the execution of 
CLOSEST-POINT, 5([0, 1]) = and L{[0, 1]) = 1 so the base case is clearly true. Suppose an interval 
Ji is split into Jn and Jj2- Because no new squeezes have occurred, 5(Jji) = S{Ji2) = S{Ji) 
and L(Jii) = L{Ji2) = L(Ji)/2. So 2^(-^')L(Ji) = 2^(-^»i)L(Ja) + 2^(-^«)L(Ji2) and the sum is not 
changed. If the interval Jj is squeezed into Jji, then S'(Jii) = 1 -|- S{Ji) and L^Jn) = L(Jj)/2 so 
2^(-^^)L(Ji) = 2^(-^»i)L(Jii) and the sum is not changed in this case either. □ 

Proposition 10. There are 0(|P| log(e~^/|P|)) samples in Phase 1. 

Proof: As noted above, we only need to count the squeezes. By Proposition IHl at the end of 
Phase 1, if Ji, . . . , Jfc contain points of P, then "Y^^^i 2^^'^^'^ L{Ji) = 1. But because no interval of 
length e or less ever appears, L(Jj) > e so X^^Li 2*^^"^'^ < e"'^ and k < \P\. Using Proposition d 
we get X]i=i S{Ji) < klog{e~^ /k) = 0{\P\ log(e~^/|P|)). Every squeeze increases S{Ji) by 1 

and no operation ever decreases it, so the number of squeezes is at most 0(|-P| log(e~^/|P|)). □ 
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Figure 4: An example C for the adaptive lower bound for n = 24 and k = 4. 



3.8 Phase 2 

If in Phase 2 a point x is sampled for which ||C(x)|| < dmin + the algorithm stops. Because 
we are giving an upper bound on the running time, we may assume that every point sampled is 
farther than dmin + £/2 from the origin. 

If closest-possible(rEj, > d^nim then will never be chosen for subdivision. This 

is because an interval around a point that is dmin away from the origin has its ellipse at distance 
at most drain from the origin and will be chosen over Thus, let us call an interval 

[xjjXj+i] alive if closest-possible(xj, Xj+i) < dmin and call it dead otherwise. No dead interval is 
ever subdivided. 

Proposition 11. If e < dmin, then when the CLOSEST-POINT enters Phase 2, there are 0{\P\) 
alive intervals. 

Proof: Let the alive intervals at the start of Phase 2 be [a^i, yi], [x2, 2/2]) • • • i [2;^, y^]. From the 
assumption above, ||C(3;j)|| > dmin + e/2 and ||C(7/i)|| > dmin + e/2. Because the intervals are 
alive, closest-possible(xi, yj) < dmin- This means that we can apply the Ellipse Lemma(2) to 
and to get closest-possible(a:;i, yj+i) < d^am — e/2. By Proposition |H1 P has a point in 

Because at most two segments of the form can overlap, and each one has at 

least one point of P, there must be at most 2\P\ of these segments. □ 

Now suppose we subdivide an interval [xi,X2] into and [2;,X2]. Because the algorithm 

is in Phase 2, closest-possible(2;i, X2) > dmin — e/2. By our assumption above, ||C(x)|| > dmin + 
e/2. Applying the Ellipse Lemma(3), we get that either closest-possible(xi, x) > dmin or closest- 
possible(a;, X2) > dmm (or both). This implies that when the interval is subdivided, at most one of 
the resulting intervals can be alive. 

Proposition 12. If e < dmin; then closest-point performs at most 0{\P\ log(e~^/|P|)) samples 
in Phase 2. 

Proof: Let /i , /2 . . . , /fc be the lengths of the alive intervals at time t. Define p{t) = Yli=i ^oS2 (2^i /^) • 
Because no interval of length 2e or less is ever subdivided, li> e and so each term in the sum is at 
least 1. At every subdivision, an alive interval is replaced with at most one alive interval of half the 
length; therefore, each subdivision decreases p{t) by at least 1. This implies that the total number 
of subdivisions cannot be greater than p(io) where to is the time when the algorithm enters Phase 2. 
Now consider the situation at time to- The total length of the alive intervals is at most 1, so we 
have X]i=i 2^i/e < 2e~^ . Applying Proposition [3 to this inequality and to the definition of p(to), we 
get p(to) < k\og{2e-^ /k). By Proposition CH k = 0{\P\), so we get p(to) = 0(|P| log(e-Vl-P|)), 
which means there are at most that many samples of C in phase 2. □ 
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3.9 Analysis Conclusion 



Theorem 1. // on a problem instance with e < (imin. we let n = OPT(e) log(e~"'^/OPT(e)), algo- 
rithm CLOSEST-POINT uses 0{n) samples and 0(n log n) additional time, where the constant in the 
O notation is independent of the instance. 

Proof: Combining Propositions IIUI 1121 and IHl we get that the number of samples the algorithm 
makes is 0{n). Because the samples are stored in a priority queue, which may be implemented as 
a heap, it takes O(logn) time to insert or extract a sample. Hence the algorithm uses 0(n log n) 
time for the heap operations. □ 

This theorem does not hold for d^\^ < e because then the condition a < d would not be satisfied 
when we invoke the Ellipse Lemma. The best conclusion we can make about the running time of 
CLOSEST-POINT when dmin < £ is that the number of samples is 0(l/e). Below we prove that it is 
impossible to do better with respect to OPT. 

4 Lower Bounds 

4.1 Worst-Case Lower Bound 

As mentioned in the introduction, in the worst case, we cannot do better than the trivial algorithm: 

Theorem 2. For any e > 0, there is a problem instance of nearest-point- on- curve on which any 
algorithm requires samples. 

Proof: Suppose we are given e. Let C be the constant "curve", C{x) = p for all x £ [0,1] with 
\\p\\ > e. Now, for any interval [xi,X2] ^ [0,1], Ec{xi,X2) is a circle centered at p whose radius 
is (x2 — xi)/2. This means that in any proof set, every two points are less than 2e~^ apart in the 
parameter space, so the OPT for this problem is G(e~^). □ 

4.2 Adaptive Lower Bound 

We prove that closest-point is optimal with respect to the number of samples of C. 

Theorem 3. For any algorithm, and for any k £ N, and any e £ {0,1/k), there is a problem 
instance with OPT = 0{k) on which that algorithm requires Q{klog{e~^ /k)) samples. 

Proof: Let e and k be given. We will construct a problem instance family for which k = r2(0PT) 
and the number of samples required by any algorithm is i^{klog{e~^ /k)) on at least one instance 
of that family. 

Let n = Divide the parameter space into n equal regions and group them into k groups 

of n/k regions each. In each group, let the curve have one spike in some region (and be flat in the 
other regions of that group). Let fe — 1 of the spikes point up, and let the remaining spike point 
down. See Figure |3 The origin is far below the curve and e is less than the height of a spike, so 
that the only solutions to a nearest-point-on-curve instance of this form are on the spike pointing 
down. Because an omniscient adversary may force the last spike the algorithm examines to be the 
one pointing down, and the algorithm cannot determine whether a spike points up or down without 
sampling on it, the algorithm must find every spike. Note that if x is a point in parametric space 
that corresponds to the boundary between groups, C{x) does not depend on where the spikes are 
chosen. Moreover, sampling inside one of the k groups (and not on a spike) only gives information 
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about whether the spike in that group is to the left or to the right of the point sampled. This 
implies that the algorithm must perform a binary search on each of the k groups. The minimum 
number of samples to do this is indeed Q{klog{n/k)). 

To show that k = Q{OPT), we note that because the curve is piecewise-linear (and parametrized 
at unit speed), placing a point at every corner gives a proof set for any e, because that completely 
determines the curve. Each spike has 3 corners and there are possibly two more endpoints, so 
OPT <3k + 2. □ 



4.3 Lower Bound for dmm < £ 

Theorem 4. For any algorithm and for any e > 0, there is a problem instance with dmin ^ e such 
that the algorithm requires Q{OPT{e)e~^) samples to solve it. 

Proof: Because dmin < e, OPT(e) = 1. To define C, let us split [0, 1] into e~^/4 intervals of width 
4e. Fix one of these intervals, J = {xi, xi + 4e). For x J, let C(x) = (0, 2.5e). For x G J let 

„, . _ J (0, 2.5e — (x — xi)) for x < xi + 2e 

~ \ (0, 2.5e - (xi + 4e - x)) for x > xi + 2e. 

Informally, one of the intervals has a spike of height 2e pointing at the origin. Now, at x = xi + 2e, 
C(x) = (0,e/2), so dniin = and OPT = 1. The only valid outputs on such a problem instance 
are points on the spike. Because sampling the curve anywhere except J gives no information on 
the location of the spike, which could be in any of e~^/4 possible intervals, an algorithm is forced 
to do a linear search that requires Q{e~^) samples. □ 

These lower bounds also work for randomized algorithms, because the reductions are from 
linear search and binary search, problems for which randomized algorithms can do no better than 
deterministic algorithms (up to constant factors). 



5 Farthest-Point-on-Curve 

It is natural to consider the symmetric problem of finding a point on C whose distance to a 
given point is within e of the largest possible. It is straightforward to modify closest-point to 
FARTHEST-POINT, which solves the farthest-point-on-curve problem. The first two lower bounds for 
nearest-point-on-curve hold for farthest-point-on-curve as well. The analysis is also easy to carry 
over to FARTHEST-POINT, with one exception: the natural "inversion" of the Ellipse Lemma is false. 
Figure El illustrates this. Nevertheless, the algorithm running time is the same (to within a constant 
factor) because we can prove a modified inverted Ellipse Lemma. Note that farthest-possible(xi, X2) 
(the analogue of closest-possible(xi, X2)) refers to the maximum distance from a point in Ec{xi,X2) 
to the origin. To simplify the proof, we impose an extra condition that ||C(xj)|| < d — a, which was 
required only for i = 2 in the original lemma. 

Inverted Ellipse Lemma. Let < xi < X2 < X3 < X4 < 1. Also, let d, a G M such that 

< a < d. // farthest-possible(xi, X2) > d, farthest-possible(x3, X4) > d, and ||C(xi)|| < d — a for 

1 G {1,2,3,4}, then farthest-possible(xi, X4) > d + |a. 

First, we state and prove what is essentially a special case of this lemma and then, to prove the 
lemma, we will reduce the general case to this special case. 
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Figure 5: A counterexample to the inverted Ellipse Lemma (ellipses are to scale) 



Proposition 13. Let a > 0, let A and B he points such that \\A\\ = \\B\\ = 1 and let Q be the 
intersection of the bisector of angle AOB with the circle ||7;|| = 1 + a. Let P be the intersection of 
ray OB with circle \\v\\ = 1 + 8a/5. Then \\A - Q\\ + \\Q - B\\ > \\A - P\\. 

Proof: See Figure IHl First, note that \\A — Q\\ = \\Q — B\\ so we need to show that 2\\A — Q\\ > 
\\A — P\\. Let c be the cosine of angle BOQ. Then the cosine of the angle BOA is 2c^ — 1 by the 
double angle formula. Using the Law of Cosines, we write: 



1^ - Qll = Vl + (l + a)2- 2(1 +a)c 



\A-P\\ = Jl+{l + ^] -2{l + ^){2c^-l). 




5 J 

We prove the Proposition by expanding (2||A — — ||^ — P|p and showing that it is positive: 

8a\ ^ „ ,„ 9 , s 8a 



A\\A - - P - ^11 = 4 + 4(1 + ay - 8(1 + a)c - 1 - i^l + — j + 2(2c^ _ l) j^l + _ 
which simplifies to 

4(c-l)' + ^ + y(c-l)(4c-l) (1) 

Notice that is quadratic in c and the term has a positive coefficient. This means that it 
has a single global minimum as c varies. Because we are trying to prove that is positive, it is 
sufficient to show that the expression takes on a positive value when c is at the minimum. To find 
the minimum, we differentiate with respect to c to get 

64a + 40 

c- 8a - 8, 
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which is equal to for c = f^ff • Substituting it back into (Q), we get 

36a2 36a2 240^(15 + 12a) _ 2880^ 
(5 + 8a)2 ^ ^ 5(5 + 8a)2 ~ 25(5 + 8a) " ' 

as desired. □ 

Proof of Inverted Ellipse Lemma: As in the original Ellipse Lemma, assume without loss of general- 
ity that X2 = X3. Let us also assume that the condition on the points is tight (it is straightforward 
to reduce to this case using the triangle inequality), that is, ||C(xi)|| = ||C(x2)|| = ||C(x4)|| = d — a. 
Once again, let P be the intersection of the circle \\v\\ = d + |a and the ray from the ori- 
gin through (7(3:2). We need to show that \x2 — xi\ > ||C(xi) — P|| and the symmetric case 
(|x4 — X2I > ||C(x4) — P\\) will follow. From combining them, we will be able to conclude the 
lemma. When applying Proposition 1131 by scaling the entire picture, we can assume that d — a = 1 
and so d = 1 + a. 

In general, we only need to consider the case when the ellipse Ec{xi, X2) is as small as possible, 
that is, it is tangent to the circle \\v\\ = d. Let Q be a point of tangency. Then (from the 
optical properties of ellipses) the origin O is on the bisector of angle C(xi), Q, C(x2). Therefore, 
if we reflect the line Q,C{x2) off the line QO^ C{xi) must be on an intersection of the reflected 
line and the circle ||f || = d — a. Let D be the other intersection. Because the reflection of 
C(x2) is both on the circle ||f || = d — a and on the reflected line, it could be either C(xi) or 
D. If it is C(xi), we can use Proposition [T3l with A <— C(xi) and B ^ C(x2) to conclude that 
||C(xi) -Qll + ||Q-C(x2)|| > ||C(xi) Otherwise, is the reflection of C(x2) across QO and 
there are two possibilities: 

1. If D is between C(xi) and Q, then we apply Proposition El with A <^ D and B <— C(x2) to 
get \\D - Q\\ + \\Q - C(x2)|| > \\D - P\\. We add \\C{xi)-D\\ to both sides and apply the triangle 
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Figure 7: Illustration for the second case of the proof 



inequality to get 

\X2 - xi\ > \\C{xi) - Q\\ + \\Q - C{x2)\\ = \\C{xi) - D\\ + \\D - Q\\ + \\Q - C{x2)\\ > 
> \\C{xi) - D\\ + \\D - P\\ > \\C{xi) - P\\. 

2. Otherwise, C(xi) is between D and Q. See Figure [3 Let Q' be the reflection of Q across the 
bisector of angle C(xi),0,C(x2). Then, \\C{xi)-Q'\\ = ||C7(2;2)-Q|| and \\C{x2)-Q'\\ = \\C{xi)- 
Q\\. Now the line segment C{xi), Q' has a second intersection with the circle \\v\\ = d — a (besides 
the endpoint C{xi)) because this line segment is the result of reflecting line segment DQ (which 
contains C{xi), by assumption) first across OQ and then across the bisector of C(2;i), O, C(x2). We 
will call this intersection F. F is the image of C{xi) under the two reflections described above, which 
implies that ||C(xi) — Q\\ = \\F — Q'\\. This means that ||F — Q'\\ = \\Q' — C(x2)|| so we can apply 
Proposition ini with F, B ^ C(x2), and Q ^ Q' to get \\F-Q'\\ + ||Q'-C(x2)|| > \\F - P\\. 
We add ||C(xi) — F|| to both sides to get 

\X2 - xi\ > \\C{xi) - Q'W + HQ' - C{x2)\\ = \\C{xi) - F\\ + ||F - Q'\\ + \\Q' - C{x2)\\ > 
> ||C(xi)-F|| + ||F-P|| > \\C{xi)-P\\, 

as desired. □ 

Because the Inverted Ellipse Lemma has a weaker conclusion, in terms of the constant, than the 
original Ellipse Lemma, the analogue of Proposition |H1 based on the Inverted Ellipse Lemma states 
that OPT(5e/8) < 20PT(e). This means that in order to get that OPT(3e/8) = 0(OPT(e)), 
which we need for the analysis of Phase 2, we need to apply the analogue of Proposition IHl three 
times (because (5/8)^ < 3/8). 

The analysis of farthest-point-on-curve does not have the problem that nearest-point-on-curve 
has when (imin < £• Every time the Inverted Ellipse Lemma is used in the transformed proof, the 
condition that a > d holds regardless of the curve or e, unlike in nearest-point-on-curve. 
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Theorem 5. On farthest-point- on- curve problem instance {C,e), letn = OPT(e) log(e-VOPT(e)). 
Then algorithm FARTHEST-POINT uses 0{n) samples and 0(n log n) additional time. 

6 Relative Error 

We now examine modifications to our problems in which the goal is to guarantee a relative error 
bound instead of an absolute error bound. Specifically, for the nearest-point-on-curve problem, the 
objective is a parameter x such that ||C(x)|| < (1 + e)dmin; and for farthest-point-on-curve, we 
need ||C(2;)|| > (imax/(l + We require that a nearest-point (farthest-point) problem instance has 
dmin (dmax) nonzero, because otherwise the problem is unsolvable. It turns out that simple modi- 
fications to the absolute-error algorithms analyzed above yield adaptive relative-error algorithms. 
For proving an upper bound on the number of samples used by the algorithms, we focus on the 
nearest-point problem; for farthest-point, the upper bound analysis is analogous. 

We start by defining a proof set for a relative-error nearest-point problem instance. Let P be 
a set of samples of C that includes and 1, let Up he the distance from the nearest point of P 
to the origin, and Lp be the distance from the nearest ellipse around adjacent points of P to the 
origin. We say that P is a proof set for the relative-error problem instance (C, e) if Lp > and 
Up/Lp < 1 + e. It is easy to show the analogue to Proposition 0J that a proof set for relative error 
is required for a relative-error algorithm to guarantee correctness. So relative-error OPT is the size 
of a smallest proof set, minus at most 2 to account for the fact that including and 1 may not be 
necessary. 

To modify the absolute-error algorithm closest-point, first note that as it executes, dmin is 
an upper bound on dmim and the top element of Q is a lower bound on dmin- Let us call these 
values U and L, respectively. The termination condition in line 5 is that U — L < e. If we replace it 
by the condition that L > (to prevent division by zero) and U/L < 1 -\- e, we get a relative-error 
algorithm. 

Theorem 6. The modified algorithm for the relative- error nearest-point-on-curve problem uses 
0(OPT • log(2 + (1 + £-1) • d^^^yOPT)) samples. 

Proof: Let Eabs = ^t+t'- Notice that ii U — L < Eabs, then because L < dmin < U, 



U U dmin 

— < < — = — = 1 + e. 

L U — Sabs dmin ~ £abs 1 ~ 14:^ 

So the relative-error algorithm with error e terminates no later than a hypothetical execution of the 
absolute-error algorithm would with error Eabs- By Theorem Q we know that such an absolute-error 
algorithm terminates after at most 0(OPTabs ■ log(2 + eAils/OPTABs)) samples, where OPTabs is 
the absolute-error OPT for £abs- We now have an upper bound on the running time of the modified 
relative-error algorithm in terms of OPTabs- To complete the proof, we need to show a lower 
bound on OPT in terms of OPTabs- 

In a relative-error proof set P, Lp > dmin — ~^f7^> because otherwise, 

Up ^ dmin dmin 1 .. 

C(mm jq:^ ^ l+e 

So if we take a proof set for relative error e and add a sample at distance dmin from the origin, we 
obtain a proof set for absolute error Eabu- This proves that OPT(e) + 1 > OPTabs- On the other 
hand, if we have an absolute-error proof set P for Sabs; we have Up — Lp < Eabs, so Up/Lp < (l+e). 
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Figure 8: An example curve segment on which the proof of Theorem |7| is based 



which imphes that it is also a relative-error proof set for e, and so OPT(e) < OPTabs- Therefore, 
the relative-error algorithm performs O ^OPT(e) • log ^2 + -^^^-^ j OPT(e)^^ samples. □ 

We modify the construction used in proving Theorem |31 to prove a lower bound for the relative- 
error problem. 

Theorem 7. For any algorithm and for any < e < 1 and A; G N, there is a problem instance 
with OPT = 0{k) on which that algorithm requires ^l{klog{e^^)) samples to solve the relative- error 
problem. 

Proof: Consider a piecewise-linear curve segment as shown in Figure |HJ Because such a segment is 
piecewise-linear, 5 samples are sufficient to obtain all information about it. We would like to show 
that for some combinations of S, L, and D, the only solutions to the relative-error problem are on 
the spike and it takes logarithmic time to find it. 

In order for the only solutions to the relative-error nearest-point problem to be on the spike, 
the distance from O to the tip of the spike has to be smaller than D/(l + e). The distance from 
the tip of the spike to O is maximized when the spike is at one of the endpoints of the curve 
segment. In this case, the distance from the tip of the spike to O is + — S. So we need 

D > {1 + e)(/D2 + l74 - S), which is equivalent to S/L > y^iD/Lf + l/A - D/{L + Le). If 
D/L = — , ^ , the inequality becomes: 

S I 1 1 1 Ve2 j_ 2e 

L ^ V 4e2 + 8e ^ 4 " (1 + e){2Ve'^ + 2e) ~ 2e + 2 ' 

So if we choose S = L-^e, the above inequality is satisfied (because {2e + 2)y/e = V 4e^ + Se^ -I- 4e > 

V^^T2^). 

Therefore, we can construct a curve segment of arbitrarily small length L + 2S with the only 
solutions to the nearest-point problem on a spike of size 25*, which is no more than 2L\fe. Sampling 
on the curve segment but not on the spike only gives information whether the spike is to the left or to 
the right of the point sampled. Therefore, a binary search taking 0(log((L + 25)/5)) = r2(log(l/e)) 
steps is necessary to find the spike. 

To construct the curve, simply paste k copies of curve segments, as described above (they may 
overlap), except make — 1 of the spikes point away from the origin and only one point toward it. 
Because the length of each curve segment can be arbitrarily small, the total length can be made 
exactly 1 (and therefore, appropriately parameterized, is a valid input). The only solutions are on 
the spike pointing toward the origin. As in the argument for Theorem|21 a binary search is required 
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to find each spike and a linear search on the curve segments is required to find the spike pointing 
toward the origin, giving a lower bound of r2(A; log(l/e)). On the other hand, OPT < 5k. 

For farthest-point, the construction is analogous to the one above, but "flipped". On the curve 
segment on which the solution is located, the spike points away from the origin. To ensure that 
the only solutions are on the spike, the distance from the tip of the spike to O has to be at least 
+ (L/2)2. This distance is minimized when the spike is in the middle of the curve segment 
and the distance from the tip to O is D + S*. Thus, we need D + S > (1 + e)-sjD'^ + which is 

the same as SjL > (1 + e) {D/LY + 1/4 — D/L. Notice that the right hand side is simply (1 + e) 
times the right hand side of the analogous inequality for nearest-point. Therefore, \i D /L is as for 
nearest point and S = L{1 + e)\/e (which is still 0{Ly/e)), the inequality is satisfied and a binary 
search on each curve segment requires r2(log(l/e)) samples. □ 

The upper and lower bounds for the relative-error problem do not match. We leave open the 
problem of finding an optimally adaptive algorithm in this setting. 

7 Conclusion 

The results in this paper give asymptotically tight bounds on the absolute-error nearest-point-on- 
curve and farthest-point-on-curve problems in the adaptive framework. We also show almost tight 
bounds in the relative-error setting. We believe that a similar analysis can provide insight into 
the adaptive performance of algorithms for other curve problems based on Proposition ^ including 
those described in |GW90j . We plan to carry out this analysis in the future. A more difficult open 
problem is generalizing Proposition ^ from one-dimensional curves to two-dimensional surfaces in 
a way that allows algorithms based on the generalization. 
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